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REFLECTION  OF  A  SHORT  NARROW  LIGHT  PULSE 
FROM  A  SCATTERING  AND  ABSORBING  OCEAN* 

Vladimir  I.  Haltrin 

Naval  Research  Laboratory,  Ocean  Sciences  Branch,  Code  7331,  Stennis  Space  Center,  MS  39529-5004,  USA 
Phone:  601-688-4528,  fax:  601-688-5379,  e-mail:  <haltrin@nrlssc.navy.mil> 

ABSTRACT 

A  very  simple  analytical  expression  for  the  shape  of  a  pulse  reflected  from  scattering  and 
absorbing  seawater  is  obtained.  The  resulting  equation  can  be  used  for  algorithms  connected  with 
the  rapid  assessment  of  optical  water  properties  from  remote  platforms. 


1.0  INTRODUCTION 

Even  though  there  is  a  multitude  of  publications  devoted  to  the  propagation  of  a  light  pulse  through  a 
scattering  and  absorbing  medium  like  seawater  (see  Refs,  in  Dolin  and  Levin,  1991),  very  few  of  them  are  practical 
enough  to  be  used  in  real-time  detection  algorithms,  mainly  because  of  the  complexity  of  the  resulting  equations. 
Instead,  inadequate  and  oversimplified  expressions  are  coded  in  many  real  detection  programs.  In  this  article  an 
attempt  is  made  to  derive  a  very  simple  analytical  expression  for  the  spacial-temporal  shape  of  a  light  pulse 
propagating  in  seawater.  The  resulting  equation  has  a  very  simple  form  and  depends  parametrically  on  the 
characteristics  of  the  emitter  and  detector,  as  well  as  the  inherent  optical  properties  of  water. 

1.1  Emitter  and  Detector  Parameters 

First,  let  us  specify  the  parameters  of  an  emitter  and  detector.  For  simplicity,  a  Lambert-Gaussian  detector  is 
assumed.  Such  a  detector  adequately  emulates  the  majority  of  real  detectors.  The  sensitivity  of  this  detector  is 
regarded  as  Lambertian,  Le„  it  is  independent  of  the  angle  of  incidence  of  light.  The  sensitivity  of  the  detector 
surface  declines  with  the  distance  p  from  the  center  of  the  detector  according  to  the  Gaussian  law: 


4p^ 


Ttp 


(1) 


where  is  the  sensitivity  radius. which  is  defined  by 

Po=^]up)p"dp,  (2) 

0 

and  is  the  detection  efficiency  of  the  receiver 

ko  =  2KjTo(p)pdp.  (3) 

0 

We  assume  an  emitter  that  is  Gaussian  both  over  the  angle  and  over  the  distance  from  the  center.  It  emits  an 
infinitely  short  pulse  represented  by  a  temporal  delta  function.  Such  an  assumption  is  mathematically  convenient 
because  the  response  from  any  arbitrary-shaped  pulse  can  be  calculated  by  mere  convolution  of  the  delta-shaped  pulse 
response  with  the  shape  function  of  the  real  pulse  (Morse  and  Feshbach,  1953).  The  energy  density  of  the  light  pulse 
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A,(r,n)  =  -^  +  jLp2(r,n)/7+(-nn')rfn',  (21) 

Aj(r,n)  =  ^  J Lp2 (r. n) p_ (-nn') dn'  J jt'pi(r.n);7+(-nn')r/n'.  (22) 

a,  n. 

2.2  Small-Angle  Approximation 

Let  us  choose  the  phase  function  components  and  p_  in  such  a  manner  that  their  tails  in  the  backward 
hemisphere  are  exponentially  small,  so  that 

Kl«|6i|  and  |A2|«|i22|-  (23) 

The  resulting  phase  function  from  Eqn;(10)  with  its  parU  satisfying  the  inequalities  above  still  gives  us  a 
very  satisfactory  approximation  for  a  realistic  ocean  phase  function  (Haltrin,  1984). 

Let  us  solve  Eqns.  (18)  in  the  small-angle  approximation  (Wells,  1982;  Walker,  1987;  Arnush,  1972;  Dolin 
and  Levin,  1991)  with  the  phase  functions  described  above.  We  should  also  make  the  following  simplifications  that 
are  typical  for  the  small-angle  approximation: 

n  =  n^+s,  nj=(0,  0,  l-s^/2),  nV^  s  +  ,  p±(nn0->  p±(|s'|) , 

nn'  =  l-(n-n')^/2,  n-n'  =  s',  n' =  n  +  (n'-n)  =  l  +  (s -s'), 

L^,(r,n)-4L^,(z,p,s),  tp,.(r,n')^  Lp,(z,p,s-s').  . 


Now  we  have  the  approximate  system  of  equations  for  the  Laplace  transforms  of  radiances: 


Lp,(z,p,s)  =  (2,(z,p,s), 

> 

Lp2(z,p,s)  =  !22(z,p,s). 

(25) 


The  Fourier  transform  of  radiance  in  the  plane  that  is  orthogonal  to  the  direction  of  the  pulse  propagation  is: 

L^iz,p,s)  =  jjdkdqF^,UXq)e-'^^-‘'^,  (26) 

1 


/V/(z-k,q)  =  ^  JJ dpds  L,,(z,p,s)e'’^«''» 


(27) 


Now  we  have  the  following  system  of  equations  for  the  Laplace-Fourier  transforms  of  the  forward  and  backward 
radiances: 


(28) 


r  r.  1 2V  ,  ^ 

1,  2  Jdz  dq 

(z,  k,  q)  =  V,  (q)  Fp,  (z,  k,  q) -1- y_  (q)  Fp2  (z- k,  q) , 

(  12'!^  1  ^ 

[y  2  Jdz  J 

F^2(z,k,q)  =  n(q)F^2(z.k,q)+y.(q)F^,(z,k,q), 

where 


Jjomp±i0)de , 


fc(i-/) 


4;r 


(29) 
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and  /o  is  the  zero-order  Bessel  function.  Retaining  only  terms  proportional  to  q^,  we  get; 


K(q)  =  fe(l-/)-fe(e^>q'/4,1 
V_(q)  =  bf-b{(7t-ef}q^/4,\ 


where  the  angular  brackets  (...)  denote  averaging  over  the  phase  function  given  by  Eqn.  (10)  according  to  the  rule: 

1  r 


Equations  (28)  for  the  light  pulse  components  should  satisfy  the  following  boundary  condition: 

=  Fpo(k,q)  =  Po  expf 


n 


(31) 


(32) 


Now  we  can  estimate  the  terms  in  the  left  part  of  Eqn.  (30)  which  are  proportional  to  s^: 

is^|L,~(0^)cL,.  (33) 

At  the  same  time  our  corrections  due  to  the  phase  function  in  Eqn  (27)  have  the  following  order  of  magnitude: 


(34) 

So,  if  the  condition  (P^)«  holds,  all  terms  which  are  proportional  to  in  the  left  part  of  Eqn.  (28)  may  be 
neglected.  With  this  in  mind,  Eqns.  (24)  acquire  the  following  form: 


+  AqO  ^pi(z.k,q)-  {bf-P^q^)  p2(z,k,q)  =  0, 

dz  dq  ) 

-(^/-Aq')  Pp,(z.k,q)+  Pp2(z.k»q)=0- 


(35) 


where  a-a  +  bf  piv,  >Ia  and  ^2  boundary  condition  for  Eqns.  (35)  is 

given  by  Eqn.  (32). 

Next,  let  us  represent  the  downward  pulse  radiance  as  a  sum  of  the  unscattered  part  (the  source)  and  the 
scattered  part  : 

P^,(z,k,q)  =  P«(z,k,q)+P;,(z,k,q).  (36) 


The  backward  pulse  radiance  consists  only  of  scattered  radiation  The  unscattered  forward 

pulse  radiance  F^^  satisfies  the  following  propagation  equation: 


^  d  ,  d  . 
dz  oq  j 


Pp^(z, k,q)  =  0. 


(37) 
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with  the  same  condition  on  the  boundary  given  by  Eqn.  (29).  The  solution  to  Eqn.  (37)  is  given  by  the  expression: 


F;?(z,k,q)  =  F;(k,q  +  fe)exp  -oz - ^  J(q k7]f  dr\  . 


This  expression  can  also  be  represented  by  the  following  analytical  formula: 


F;?(z,k,q)  =  Foexp[-^-^^-fa-h-^kq  +  Aq')2-f:^k'-h)3,kq]z^-A^z']  •  (39) 


The  scattered  parts  of  the  forward  and  backward  radiances  of  the  pulse  are  satisfied  by  the  following  equations: 


-•^+k^+«  +  l3,q^j  F;2(z,k,q)  F^^(z,k,q), 

•^-k^  +  a  +  Aq^j  F/i(z,k,q)  , 


with  the  boundary  conditions 


F/i  (z,  k,  q)[^^  =  0 ,  lim,_^  F;;.  (z,  k,  q)  =  0 ,  /  =  1, 2  . 


The  next  step,  according  to  the  method  given  by  Both  (1929),  Snyder  (1949)  and  Romanova  (1968),  involves 
the  following  substitutions: 

q  =  g-kz,  F^.(z,k,q)  =F^.(z,k,g-kz)  =  <I>',(z,k,g),  (42) 

which  convert  Eqns.  (40)  into  the  following  one-dimensional  system  of  two  equations: 

-■ ^  +  a+A(g-kz)^  4>p2(z,k,g)  =  [Z>/-A  (g-kz)^]  F®(z,k,g-kz), 

=  f  1  r 

— -l-a+A  (g-kz)^  <f»pi(z,k,g)  =  [bf-p2  (g-kz)^]  «>p2(z>k,g). 

The  solutions  to  Eqns.  (43)  can  be  found  with  the  help  of  the  following  two  (forward,  ,  and  backward,  G_)  Green 
functions  (Vladimirov,  1971): 

G±(z,k,g)  =  //(±z)exp[TozTA}(g-kJ7)'^n  ,  (44) 


which  are  solutions  to  the  equations: 


•-t-a-i-A(g-kz)^  G±(z,k,g)  =5(z). 


where  jFf(z)  (i/  =  l,  for  2>0,  //  =  0,  for  z  <  0)  is  a  Heavyside  or  step  function  (Morse  and  Feshbach,  1953). 
The  final  solutions  to  Eqns.  (43)  are: 


(z.  k,  g)  =  J  G.  (z  - 1.  k,  g)  [bf  -  A  (g  -  k^)^  ]  F«(^,  k,  g  -  k^)  . 
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(z. g)  =  J G+  (z - k, g) [bf - ^2 (g - ] 0*2 (|, k, g)) +  C, (k, g) exp  -ocz -Pi jig- krj)^ r/?7  ,  (47) 

L  0 

where  C,(k,g)  is  an  arbitrary  function  that  is  determined  by  the  boundary  condition. 

After  simplification,  we  obtain  the  following  expressions  for  fj,, : 

Fpiiz,k,q)  =  F^(k,q  +  kz)exp[-az-  p^jiq  +  kr|p  dT^]\l  +  j  d^[bf  -  P2(q-k^f]x 

0  [  0 

1 

xexp[-A  jiq  +  kr)fdr]]  jdg[bf-p,iq-k^+kgf]&xp[2ag-pi  jiq-k7]fdr]]\, 

i-z  -  f  J 

0  z—^ 

Fp2(z,k,q)  =  F^ik,q  +  kz)e-‘‘-^  jd^[bf-p^(q-kz  +  k^)^]exp[2a^-p^jiq-kr^fdr|]  .  .  (49) 


2.3  Detector  Response  to  the  Infinitely  Short  Laser  Pulse  Reflected  from  Water 

After  the  integration  of  the  received  radiances  over  the  sensitive  area  of  the  detector,  the  relative  (normalized  by 
the  pulse  power  P^)  response  of  the  detector  placed  at  the  depth  z  will  be 


rjiiz,w)  =  - 


6iw)ve 


Tl,.4AcrJ 

IpJ 

H  \ 

Pa 

where 


2  a.  +  2  D2  (2 1'  -  w  ^)  +  Z),C7o  ^  +  2w  j 


=P^+Po+^oZ^+Avv"+-^^(4z"  +  w"),  w=yt-z. 


D,a„ 


(50) 


(51) 


772  (z,  w)  =  e(w)  exp[-a  (w  +  z)] , 


D=l  + 


■^1  +  A 
.Pa 


1+ 


F>i<y„w 
3D. 


2  J 


( 

W 

kP^j 


3*2  J 


f  \ 

IpJj 


At  z  =  0,  the  detector  response  tdthe  reflected  pulse  will  be: 


where 


1  +  flj  T  +  Cj  T 


7  =  t/to  3  fo  =(«V)‘'  ,  Go  = 


a.  = 


nD. 


8^'(pi  +  PD)[l-®oa-/)f  ’  24c^(p^+p^)[l-cOo(l-/)f  ’ 


[l-©o(l-/)][l+(p^/p^f] 

;ra)n(l-/)D, 


(52) 


(53) 

(54) 

(55) 


where /is  the  weight  coefficient  in  Eqn.  (10). 

The  values  of  r]^  calculated  according  to  Eqns.  (53)-(55)  were  compared  with  the  results  computed  with  the 
Monte  Carlo  code  by  Kattawar  (1992).  The  discrepancies  between  the  natural  logarithms  of  the  772  computed  by  both 
these  methods  for  inherent  optical  properties  and  phase  functions  by  Petzold  (1972)  do  not  exceed  15%. 
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3.0  CONCLUSIONS 

Relatively  simple  analytical  equations  derived  in  a  small-angle  scattering  approximation  have  been  obtained  for 
the  infinitely  short  light  pulse  reflected  by  seawater.  These  equations  depend  on  the  inherent  optical  properties  of 
seawater,  as  well  as  the  parameters  of  the  pulse  and  receiver.  They  can  be  transformed  into  equations  for  an 
arbitrarily-shaped  pulse  by  a  simple  convolution  procedure.  The  logarithmic  precision  of  these  equations  is  estimated 
to  be  in  tlie  range  of  15%.  The  results  of  this  paper  can  be  used  for  algorithms  connected  with  the  rapid  assessment 
of  water  properties  from  remote  airborne  and  shipborne  platforms. 
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